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1 Introduction 

The study of reaction-diffusion partial differential equations (PDE's) of the form 
du 

— (u,t) = F(u(u,t),t) + DAu(u,t), (1) 

where A denotes a diffusion operator, is central to many biological applications, 
in fields ranging from pattern formation in development to ecology. One of the 
central topics of research in this context is the question of how the stability of 
solutions of the PDE relates to stability of solutions of the underlying ordinary 
differential equation (ODE) = F(x(t),t). This paper shows that when solu- 

tions of this ODE have a certain contraction property, namely fj,2 p(Jf(u, t)) < 
uniformly on u and t, where ^2,P is a logarithmic norm (matrix measure) associ- 
ated to a P-weighted L 2 norm, the associated PDE, subject to no-flux (Neumann) 
boundary conditions, enjoys a similar property, if P 2 D + DP 2 > 0. This result 
complements a similar result shown in [T] which, while allowing norms LP with 
p not necessarily equal to 2, had the restriction that it only applied to diagonal 
matrices P. Here, P is allowed to be an arbitrary positive definite symmetric 
matrix. The paper also discusses an example of biological interest, as well as 
examples that illustrate when the results in ppj apply but the current result does 
not. 



2 Main Result 



In this section, we study the reaction diffusion PDE dU, subject to a Neumann 
boundary condition: 

Vu i -n(^,t) = o vc e an, v*e[o,oo). (2) 

Assumption 1. In ([!]) — ([2]) we assume: 
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Q is a bounded domain in R m with smooth boundary dfl and outward normal 



n. 



• F: V x [0,oo) — > W 1 is a (globally) Lipschitz and twice continuously differ- 
entiable vector field with respect to x, and continuous with repect to t, with 
components Fj/ 

F(x,t) = (Fi(x, i), • • ■ ,F n (x,t)) T 
for some functions i^: V x [0, 00) — > R, where V is a convex subset ofW 1 . 

• D = diag(di, • • • ,d n ), with di > 0, is called the diffusion matrix. 
Definition 1. By a solution of the PDE 

du 

— (u>, t) = F(u(u, t),t) + DAu(uj, t), 

Vui ■ n(f , t) = V£ G dU, Vt G [0, 00). 

on an interval [0, T), where < T < 00, we mean a function u = (u±, ■ ■ ■ ,u n ) T , 
with u: 0, X [0, T) — > V, such that: 

1. for each u) G Ci, •) is continuously differentiate; 

2. for each t G [0,T), u(-,t) is in Y, where 

Y = L-.n^v l « = ,«„), Wi GCg(n), ^(O = o, v^ean ViJ, 

where Cj| (O) is i/ie set 0/ twice continuously differ entiable functions — > 
R. 

5. /or eac/i we!!, and eac/i i G [0, T), u satisfies the above PDE. 

Theorems on existence and uniqueness for PDE's such as (pQ) — ([2]) can be found 
in standard references, e.g. [HE]. 

For any invertible matrix Q, and any 1 < p < 00, and continuous u: CI — > R n , we 
also denote ||tt||p,Q = ||Qu|| p , where (Qu)(uj) = Qu{uj) and || • || p now indicates 
the norm in LP(fi,K n ). 

Definition 2. Let (X, \\ ■ \\x) be a finite dimensional normed vector space overM 
or C. The space C(X,X) of linear transformations A: X — > X is also a normed 
vector space with the induced operator norm 

||.A||x->jc = sup 
IMIx=i 

The logarithmic norm Hx(') induced by \\ ■ \\x is defined as the directional deriva- 
tive of the matrix norm, that is, 

Hx(A)= liin ~ + hA\\ x ^ x -l), 
h— >o+ n 

where I is the identity operator on X. 
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In pp, we proved the following lemma: 

Lemma 1. Consider the PDE system {[TJ) — d2J) ^ and suppose Assumption^ holds. 
For some 1 < p < oo, and a positive diagonal matrix Q, let 

fi : = sup n Pt Q(J F (x,t)). 
(x,t)eKx[0,oo) 

Then for any two solutions u and v of ([I]) — ([2]) , 

|K,<) - v(;t)\\ P)Q < e^|K,0) - w(-,0)[| p , o . 

Before stating the main theorem of this section, we first prove the following: 

Lemma 2. Suppose that P is a positive definite, symmetric matrix and A is an 
arbitrary matrix. 

1. If H2,p{A) = n, then QA + A T Q < 2fiQ, where Q = P 2 . 

2. If for some Q = Q T > 0, QA + A T Q < 2fiQ, then there exists P = P T > 
such that P 2 = Q and fi2,p{A) < /i. 

Proof. First suppose fj,2 t p(A) = fi. By definition of \i: 

^ (PAP~ X + (PAP- 1 )^ < ill. 

Since P is symmetric, so is P , so 

PAP" 1 + P~ 1 A T P < 2(j,I. 

Now multiplying the last inequality by P on the right and the left, we get: 

P 2 A + A T P 2 < 2fj,P 2 . 

This proves 1. Now assume that for some Q = Q T > 0, QA + A T Q < 2[iQ. 
Since Q > 0, there exists P > such that P T P = Q; moreover, because Q is 
symmetric, so is P. Hence we have: 

P 2 A + A T P 2 < 2fiP 2 . 

Multiplying the last inequality by P from right and by P~ l from left, we conclude 
2. □ 

Remark 1. Observe that for Q > 0, 
1. 

QA + A T Q < fiQ QA + A T Q < /3I, 

where j3 = fj,X and A is the smallest eigenvalue of Q. 
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2. 

QA + A T Q < pi => QA + A T Q < jQ, 

where 7 = — and A' is the largest eigenvalue of Q. 
A 

Theorem 1. Consider the reaction diffusion system (PQ) — d21) and suppose As- 
sumption [1\ holds. Let 

H := sup H2,p(Jp(x,t)), 
(x,t)eVx[0,oo) 

for a positive symmetric (not necessarily diagonal) matrix P. In addition assume 
that QD + DQ > 0, where Q = P 2 . Then for any two solutions, namely u and v, 
of dU) - ©, we have: 

IK,t)-<,t)|| 2) p<e^|K,0)-t;(-,0)||2 ) P. (3) 

Proof. By Lemma [2j 

QJ F + J T F Q < 2fiQ, (4) 

where Q = P 2 . Let u, and v be two solutions of ([I]) — ([2]) and let w := u — v. 
Define 

$0) := ~(w,Qw) = -\\Pw\\ 2 2 , 

where (x.y) := f X *y. Since u, and „ satisfy ©, we have: 
Jn 

— H = (w, Q(F(u, i) - F(«, t))) + («;, QDAw) . (5) 

Since + is a positive matrix, there exists a positive, symmetric matrix 
M such that QD + DQ = 2M T M. 

{w,{QD+DQ)Aw) = 2(w,M T MAw) = 2(Mw,AMw) = -2(V(Mw),V{Mw)). 

(6) 

The last equality holds by the Mean Value Theorem and the fact that Mw satisfies 
the Neumann boundary conditions. 

By Green's identity and using DQ = (QD) T , we have: 

(w, (QD - DQ)Aw) = {w,AQDw) - (QDw,Aw) = 0. (7) 
Equations © and ((7J imply: 

(w,QDAw) = -(V(Mw),V(Mw)) < 0. 



The techniques of the proof are similar to the proof of Theorem 1, [2] 
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Now by Mean Value Theorem for integrals, and using (J3|), we rewrite the first 
term of the right hand side of (0) as follows: 

(w,Q(F(u,t) - F(v,t))) = w T (oj,t)Q(F(u(u,t),t) - F(v(u,t),t)) doj 

Jn 

T f 

w (u,t)Q / Jp{v(uj,t) + sw(ui,t),t) • w{uj,t) ds du 
iq Jo 



< fx I ds w Qw doj 



UU 

n 

T, 



Therefore 



fj> w Qw doj 
Jn 



— H < 2n*{w) 



This last inequality implies ([3]). □ 

Corollary 1. In Theorem^ if fj, < 0, then §1$) — is contracting, meaning 
that solutions converge (exponentially) to each other, as t -> +oo in the weighted 
L>2,p norm. 

Example 1. In pQ we studied the following system: 

xt = z — 5x + kiy — /c2(Sy — y)x + d\Ax 

yt = -hy + k 2 (S Y - y)x + d 2 Ay, 

where (x(t),y(t)) £V= [0,oo) x [0, 5y] for all i > (V is convex), and Sy, k%, 
k 2 , S, d\, and d 2 are arbitrary positive constants. 

In [3], it has been shown that for p = 1, there exists a positive, diagonal matrix 
Q such that for all (x,y) 6 V, Hi } q( Jf(x, y)) < 0; and then by Lemma[3]one can 
conclude that the system is contractive. 

We showed that for any positive, diagonal matrix Q and any p > 1, there exists 
(x,y) E V such that // Pi g(Jp(x, y)) > 0, where 

F = (z - Sx + kiy - k 2 (S Y - y)x, -hy + k 2 (S Y - y)x) T , 

and 

( —5 — a b 

Jf = 



a —b 

with a = k 2 (Sy — y) € [0, k 2 Sy] and b = ki + k 2 x G [ki, oo). 

Now we show that there exists some positive, symmetric (but non-diagonal) ma- 
trix P such that for all (x,y) G V, ^ 2 ,pJp{x,y) < 0. Then for those di and d 2 
that satisfy P 2 D + DP 2 > 0, where D = diag(di, d 2 ), Corollary Q] concludes that 
the system is contractive. 

"1 1" 

} 1. 

Note that Q is symmetric and positive (because q > 1). 



Claim. Let Q 



where q > 1 + — . Then QJ F + (QJ F ) T < 0. 



5 



Proof of Claim. We first compute 



QJf 



"l 1" 




—5 — a b 




.1 ?. 




a —b 





-5 
-5+(q-l)a -b{q-l) 



So 



QJ F + ( J F Q) 5 



-25 -8 + (q-l)a 

-8 + (q-l)a -2b(q-l) ' 

To show QJf + < 0, we show that det (QJf{x, y) + Jp(x, y)Q) > for all 
G V: 



det + JJQ) = 4<»(g - 1) - (-5 + {q- l) a y . 

Note that for any q > 1, /(a) := (—5 + (g — l)a) 2 < J 2 on [0, k 2 Sy], and <?(&) := 
A5b(q — 1) > 4:5k\(q — 1) on [ki,oo]. So to have det > 0, it's enough to have 

<5 2 5 

45fci(g - 1) - <5 2 > 0, i.e. g- 1 > — — , i.e. q > 1 + — . □ 

Now by Remark [J and LemmaH for P = y/Q, /-i 2j p(Jf(x, y)) < 0, for all (x,y) G 
V. 

Example 2. We now provide an example of a class of reaction-diffusion systems 
Xt = F(x) + DAx, with x G V (V convex), such that for some positive definite, 
diagonal matrix Q, and for all x G V, ^i,q(Jf( x )) < (and hence by Lemma HI 
these systems are contractive), yet for these systems, H2,p(Jf(x)) -ft 0, even for 
any positive definite, symmetric (not necessarily diagonal) matrix P. Consider 
two variable systems of the following type 



xt = -fi(x) + gi(y) + diAx 
Vt = f2{x) - gi{y) + d 2 Ay, 



(8) 
(9) 



where d\, d 2 are positive constants and (x, y) G V = [0, oo) x [0, oo). The functions 
fi and gi take non- negative values. Systems of this form models a case where x 
decays according to /i , y decays according to g 2 , and there is a positive feedback 
from y to x (g%) and a positive feedback from x to y (f 2 ). 

Lemma 3. In system (G))-(I2)) ; let J be the Jacobian matrix of 

(-/i(z)+^fo),/ 2 (s)-^(|/)) T . 

In addition, assume that the following conditions hold for some A > 0, and \l > 
and all (x, y) G V: 



1. -fi(x) + \\&(x)\<-n<0; 



2- -g'M + j\g'M\ < -m<o ; 



3. for any p G 



lim 

y— >oo 



(g'i(y) -pog' 2 {y)Y 
g' 2 (y) 



oo. 



6 



Then 



1. for every (x,y) E V, fii ; Q{J{x,y)) < 0, where Q = diag(l,A); and 



2. for each positive definite, symmetric matrix P, there exists some {x, y) G V, 
such that n,2,p(J{x,y)) > 0. 

Proof. The proof of /j>i,q( J(x, y)) < is straightforward from the definition 
of fj,i Q and conditions 1 and 2. Now we show that for any positive matrix 

P = , there exists some (xo>2/o) G V such that \x 2 p{J{xq, yo)) > 0. 

p P2 

By Lemma [21 it's enough to show that for some (xo>2/o) £ V, PJ(xo,2/o) + 
J T (xo,yo)P -jt 0. We compute: 



PJ 
Therefore, 
PJ + (PJ) T 



Pi 


p 




p 


P2_ 





ffa) -9' 2 (y). 



-pifiix) +pf' 2 {x) pig[(y) -pg' 2 (y) 
-pf[(x) +p 2 f 2 (x) pg[(y) -P2g' 2 (y) 



2(-pif[(x) +pf 2 (x)) Pl g[(y) - pg' 2 (y) - pf[(x) + p 2 f' 2 {x) 

pig'iiv) -pg' 2 {y) - pf[(x) + p 2 f 2 (x) 2 (p9i(y) -P2g' 2 {y)) 



Now fix io G [0, oo) and let 

A := 2 (-pif[(x ) +pf 2 (x )) and B := -pf[(x ) + p 2 f 2 {x ). 

Then 

det (PJ+ (Pjf) = 2A {pg'M - P2 g' 2 {y)) - ( Pl g[(y) - pg' 2 (y) + B) 2 . (10) 
We will show that det < 0. Dividing both sides of (fTUI) by pfg' 2 (y), we get: 
det {PJ + (PJ) T ) 2A {pg[{y) - p 2 g' 2 {y)) {g[{y) - p g' 2 {y) + B'f 



Pi9 2 {y) 



plg' 2 {y) 
(g[{y) - pog' 2 {y)) 2 



g'M 



2B' 9 M + 2B' P0 



g'M 



g'M 



g'M 

B 

Pi Pi 

(Note that p\g' 2 {y) > because by condition 2, g' 2 > /i > 0, and P > implies 
Pi + 0.) 



where po = — , A' 
Pi 



2A 

— , and B 



By condition 2, < 



g'M 



< A < oo for all y. Now using condition 3, we can find 



g'M 

y large enough such that det < 0. 

Since det {PJ{xo, yo) + {PJ{xo, yo)) T ) < for some (xo,yo) £ V, the matrix 
PJ + {PJ) T has one positive eigenvalue. Therefore PJ + {PJ) T -ft 0. □ 
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As a concrete example, take the following system 

xt = —x + y 2+€ + diAx 

y t = 5x- (y 3 + y 2+e + dy) + d 2 Ay, 

where < 6 < 1, < e <C 1, d, di, and d 2 are positive constants and (x,y) £ 
V = [0,oo) x [0,oo). 

In this example we show that, the system is contractive in a weighted L 1 norm; 
while for any positive, symmetric matrix P, and some (x,y) E V, fJ<2 pJf(x,v) ft 
0. To this end, we verify the conditions of Lemma [3j 

For any (x, y) E V, we take in LemmaEl A = 1, and any [i E (0, min{d, 1 — 5}): 

1. -1 + 5 < 0, because < 8 < 1. 

2. - (3y 2 + (2 + e)y 1+e + d) + (2 + e)y 1+t = -3y 2 -d<-d<0. 

3. For any po £ R, 

((2 + e)^ - po (3y 2 + (2 + e)^ + d)) 2 

lim — ^ 7- r -r— : = CO. 

y^oo 3y 2 + (2 + e)y 1 + e + d 

So the conditions in Lemma [3] are verified. □ 

3 Diffusive interconnection of ODEs 

In this section, we derive a result analogous to that for PDE's for a network of 
identical ODE models which are diffusively interconnected. We study systems of 
ODE's as follows: 

u(t)=F(u(t))-(L®D)u(t). (11) 
Assumption 2. In (jlip . we assume: 

• For a fixed convex subset o/R™, say V, F : V N — > W lN is a function of the 
form: 

F(u) = (F(uY,... ,F(u N f) T , 

where u = [u lT ', • • • , u nT ^J , with u' 1 £ V for each i, and F : V — > R n is a 
(globally) Lipschitz function. 

• For any u £ V N we define \\u\\ q as follows: 

P 

where Q is a symmetric and positive definite matrix and 1 < p < co. 
With a slight abuse of notation, we use the same symbol for a norm in R n : 

ll^llpiQ := 1 1 1 1 - 



\n P ,c 



(\\Qu l \\ p ,--- ,\\Qu? i 
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• u: [0,oo) —7- V N is a continuously differentiable function. 

• D = diag((ii, • • • ,d n ) with di > 0, which we call the diffusion matrix. 

• L G M. NxN is a symmetric matrix and LI = 0, where 1 = (1, • • • , 1) T G R^. 
We think of L as the Laplacian of a graph that describes the interconnections 
among component subsystems. 

In pp, we proved the following lemma: 

Lemma 4. Consider the ODE system Ul\) . and suppose Assumption [H holds. 
For some 1 < p < 00, and a positive diagonal matrix Q, let 

H := sup fip } Q(J F (x,t)). 
(x,t)eVx[0,oo) 



(^p,Q i- s the logarithmic norm induce by the norm 
\\x\\ P ,Q ■= \\Qx\\ p .) 

Then for any two solutions u and v of ([lip , 



\ p< q on W n defined by 



\u(t) - v(t)\\ p , Q <e^\\u(0)-v(0)\\ p , Q . 



In this section we generalize the above result for p = 2 and any symmetric, 
positive definite (not necessarily diagonal) matrix P such that P 2 D + DP 2 > 0. 

Theorem 2. Consider the ODE system and suppose Assumption [H holds. 
Let 

fji:= SUp H2 jP (J F (x,t)), 

(x,t)ev x[o,oo) 

for a positive symmetric (not necessarily diagonal) matrix P. In addition assume 
that QD + DQ > 0, where Q = P 2 . Then for any two solutions, namely u and v, 
of (G2P, we have: 



|u(t)-u(t)|| 3 ,p<e^||u(0)-t;(0)||2 ) . 



(12) 



Before proving the theorem, we recall that if A = (ay) is an m x n matrix and 
B = (bij) is a p x q matrix, then the Kronecker product, denoted by A (g> B, is 
the mp x nq block matrix defined as follows: 



A®B :-- 



011-B 



a m \B 



where aijB denote the following p x q matrix: 



a\ n B 



Q"mnB 



aijB :- 



Qijbu . . . aijbi q 



O'ij bpi 



O^ij bpq 



The following are some properties of Kronecker product: 
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1. {A ®B){C®D) = {AC) ® (BD); 

2. {A®B) T = A T ®B T . 

3. Suppose that A and are square matrices of size n and m respectively. 
Let Ai, • • • , A n be the eigenvalues of A and fj,\, • • • , \i m be those of B (listed 
according to multiplicity). Then the eigenvalues of A <g) B are for 
i = 1, • • • , n, and j = 1, • • • , m. 

Proof of Theorem O H By Lemma [2j 

QJf + JpQ < 2/xQ, (13) 
where Q = -P 2 . Let u, and f be two solutions of (jlip and let w := u — v. Define 

:= ^w T (/at <8> Q)w. 

Note that indeed 



*m = 2 E wiT Q wi = 2 E( Pu/ ) T ( Pu/ ) = 2 E n pu; 




j = l 8=1 j=l 

Since it, and i> satisfy (llip . using the first property of Kronecker product listed 
above, we have: 

^(w) = w T (I N ® Q)(F(u, t) - F(v, t)) - w T {I N ®Q){L® D)w 
= w T (I N ® Q)(-F(u, t) - F(v, t)) - w T (L ® QD)io. 

Since Q-D + is a positive matrix, there exists a positive, symmetric matrix 
M such that QD + DQ = 2M T M. Then using 

L ® (M T M) = (I N ®M T ){L®I n )(I N ®M), 

we get: 

u; T L ^ (Q-D + DQ)iw = 2k; t (/jv ® M T ){L ® Af)to (16a) 

= 2((/^®MH T (L®/ n )((/ 7V <S>MH > 0. (16b) 

The last inequality holds because all the eigenvalues of L, and therefore all the 
eigenvalues of L<g>I n , are non-negative, by the third property of Kronecker product 
listed above. Now because L = L T , DQ = (QD) T , using the second property of 
Kronecker product listed above, we get: 

w T L ® (QD - DQ)w = w T [(L ® QD)w] - w T (L T ® {QD) T ) w (17a) 
= w T [(L ® QD)w] — w T (L ® QD) T w (17b) 
= it; T [(L ® QD)w] - [(L ® Q£>)u;] T u; = 0. (17c) 

2 The techniques of the proof are similar to Theorem 4, [2] 
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Considering the sum of Equations (|16p and (|17p . we get: 

-w T {L®QD)w < 0. 



Now by Mean Value Theorem for integrals, and using (|13f) . we rewrite the first 
term of the right hand side of (|15p as follows: 

JV 



(I N ®Q)(F(u,t)-F(v,t)) = £V Q{F{u\t)-F{v\t))% 

i=l 

S^wfQ J F {v l + sw i } t)w i ds 
i=i Jo 

fl * ,rp 

< [i \ dsj^w 1 Qi 
Jo i=i 



i=l 

,1 N 
1=1 

Hw t {In ® Q)w 



Therefore 

— H < 2 M $H. 

This last inequality implies (|i2|) . □ 
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